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A spherocylinder-like molecule with a Lennard-Jones type interaction is proposed as a model 
of smectic-j4 (Sm-A) liquid crystals, which can form a free-standing film. By means of Gibbs 
ensemble simulations, the isotropic, nematic, and Sm-A phases of the model fluid are found to 
coexist with a vapor phase; and the coexistence conditions of the liquid crystal phases with 
the vapor phase are determined. For a set of the interaction-parameters of the model molecule, 
the Sm-yl free-standing film is produced below the bulk isotropic-Sm-A phase transition tem- 
perature by using Monte Carlo simulations. The film tension of the Sm-^4 free-standing film 
is calculated and its dependencies on the temperature and on the number of molecules are 
discussed. 

KEYWORDS: liquid crystal, smectic-/4 phase, free standing film, layer-thinning transition, spherocylinder, 
Lennard-Jones potential, Gibbs ensemble simulation, Monte Carlo simulation 



1. Introduction 

Liquid crystals (LCs) in a smectic-^4 phase (Sm-A) 
form a free standing film (FSF) or freely suspended film. 
The FSF is stable below the bulk phase transition tem- 
perature between the Sm-^4 and lower ordered phase. The 
FSF consists of multiple smectic layers and thus the film 
thickness is approximately integral multiple of a single 
layer thickness. Because the FSF can be very thin up 
to the two-layer thickness, it may be strongly affected 
by its surfaces; therefore the FSF is an appropriate sys- 
tem for us to investigate the effects of free surfaces on 
an ordered phase. A surface-enhancement of Sm-^4 order 
(SESO) observed in experiments 1,2 is a phenomenon due 
to the surface effect. The SESO arises since the layer fluc- 
tuation of the Sm-^4 phase is suppressed at the free sur- 
face by the surface tension. 1-3 A so-called layer-thinning 
(LT) transition, 4-11 which occurs in fairly rare Sm-^4 ma- 
terials, is another phenomenon resulting from the strong 
surface effect. The Sm-^4 FSFs of LT materials can exist 
even at the higher temperature than the bulk isotropic- 
liquid-Sm-yl (I-A) phase transition temperature Tm 4-9 
or than the bulk ncmatic-Sm-^4 phase transition temper- 
ature T NA . 10 ' 11 The FSFs remain above T IA or T NA but 
the upper bound of the number of layers is a decreasing 
function of the temperature, so that the film thickness de- 
creases layer- by- layer as the temperature increases. Such 
layer-by-layer transitions result from a nucleation and 
growth of a dislocation loop at the innermost layer; i.e., 
the melt of the innermost Sm-^4 layer followed by the es- 
cape of LC molecules in the layer into a meniscus at the 
film border. 8-10 The stability above Tja or Tma and the 
melt of the innermost layer indicate that the Sm-yl or- 
der in LT materials is emphasized by the interfaces more 
strongly than in ordinary Sm-j4 materials. 

Although the Sm-^4 FSFs are major substances in ex- 
periments, there are no molecular models which exhibit 
the FSF. The well-known Gay-Berne (GB) model, 

which is an anisotropic molecule with a spheroid-like in- 
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termolecular interaction, has isotropic liquid(J), nematic 
(N), Sm-^4, and smectic- £? phases; 12 the GB model is 
frequently used as a model of uniaxial LC molecules. By 
using sets of model-parameters with which the GB model 
fluid exhibits vapor (V), I, and N phases, it was shown 
that the GB fluid forms a nematic FSF in equilibrium 
with the vapor phase. 13,14 The Sm-yl FSF, however, has 
not been found in GB fluids. The GB model also has a 
disadvantage that it does not allow the intuitive under- 
standing of the relation between model-parameters and 
the resulting potential since the GB potential depends on 
the parameters in a complicated way. It was shown 15, 16 
by means of the density functional theory (DFT) that the 
model fluid consisting of hard-spheroids with attractive 
anisotropic Yukawa potential forms a Sm-^4 FSF. In a re- 
stricted region of potential parameters, the model Sm-^4 
FSF exhibits the SESO, as in the experiments on ordi- 
nary Sm-^4 1,2 and on the LT materials. 9 Furthermore, 
refs. 15 and 16 indicate that the model exhibits the LT 
transition in the proximity of a V—I—N—Sm-A quadru- 
ple point. These DFT results are important to show the 
relation between the intermolecular interaction and the 
stability of the FSF. However, there always remains pos- 
sibility that a foreseen phase is stable, because only a 
limited number of phases is under consideration in the 
DFT. 

The main purpose of the present paper is to make a 
pair potential of a LC molecule such that the LC forms 
the Sm-^4 FSF and to investigate physical properties of 
the LC model fluid. The new LC model is introduced in 
§ 2. By using Gibbs-ensemble simulations, in § 3.1, the 
model fluid is shown to have /, N, and Sm-^4 phases in 
coexistence with the vapor phase. In § 3.2, it is shown 
that the model fluid forms a Sm-^4 FSF for a set of 
interaction-parameters. It is found that the film is un- 
stable above Tja, and hence the new LC model cannot 
be a model of LT materials. However, the film tension 
of the LC model shows some resemblance to that of LT 
materials. Such a resemblance is discussed in § 4. Some 
concluding remarks are also presented in § 4. 
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Fig. 1. (color-online) Two dimensional representation of sphero- 
cylinders. The equidistant surface (the distance is D/2) from 
a line segment of length L is the surface of a spherocylinder. 
Therefore the interaction between two hard spherocylinders is a 
function of d, the distance between the line segments defining 
the spherocylinders. The new interaction between the soft sphe- 
rocylinders introduced in this work is also defined as a function 
of d (see details in text). 




Fig. 2. (color-online) The Lennard-Jones spherocylinder inter- 
action V(ri2,ui,U2) for four typical configurations. The po- 
tential parameters are (si,E2>£'i) = (1,0.6,-1.4). Four pairs 
of segments shown in the figure indicate the configurations of 
molecules. 



2. Lennard-Jones Spherocylinder Model 

For the purpose to investigate the Sm-^ FSF, it is 
necessary to find first a molecular model which exhibits 
a Sm-4 phase in a bulk system. Furthermore, in order to 
understand the relation between properties of the pair 
potential and resulting phases, it is preferable that all 
the terms in the potential function of the model can be 
readily understood. 

In this paper, I make a molecular interaction based 
on a hard spherocylinder (SPC), with which the model 
fluid is known to exhibit the Sm-vl phase. 17 The SPC 
is a cylinder capped with two hemispheres; the diameter 
D of the cylinder and hemispheres and the length of the 
cylinder L determine the shape of the SPC. Because the 
surface of a SPC is the equidistant surface from the line 
segment defining the center of the cylinder (see Fig. 1), 
the pair potential between two hard SPCs is expressed 
using d, the distance between two line segments: The pair 
potential of the hard SPCs is infinite in the case that d < 
D and otherwise zero. The interaction between the new 
model molecules, named Lennard-Jones spherocylinder 
(LJ-SPC), is also a function of d and is a simple extension 
of the hard SPC model to a soft molecule: 

V(n 2 ,ui,u 2 ) = 4y o (r 12 ,Ui,U2) | (J^J - (^pj j, 

(1) 

where r\ 2 = r 2 —r\) T{ and iii indicate, respectively, the 
center of gravity and the unit vector along the symmetry 
axis of the ith molecule; r\ 2 is a unit vector T"i2/|'»"i2|- 
The function V introduces an additional orientational 
dependence: 

V (r 12 ,ili,U 2 ) = £l + £2P 2 {U1 ■ U 2 ) 

+ ^{P 2 {ri 2 -Ui) + P 2 {r 12 -u 2 )}, (2) 

where P 2 {x) = (3x 2 — l)/2 is the Legendre polynomial 
of second order. The terms in eq.(2) correspond to the 



leading order terms in a general spherical harmonic ex- 
pansion of a function of r\ 2 , Ui, and ii 2 . Similar expan- 
sion form of the anisotropic pair potential has been used 
in refs. 18-21 with a modified atomic LJ potential, and 
in refs. 15 and 16 with an anisotropic Yukawa attractive 
potential. The first term in eq.(2) gives uniform contri- 
bution and no additional orientational dependence on 
the pair potential. Thus, with e 2 — £3 — 0, the molecule 
has only the SPC-like orientational dependence resulting 
from the use of the distance d. In the case that e 2 > 0, 
the second term, which is a Maire-Saupe-type interac- 
tion, makes parallel molecules stable and it contributes 
to the nematic ordering. The third term makes, in the 
case that £3 < 0, the molecules prefer a side-by-side con- 
figuration rather than end-to-end configuration, so that 
the term induces a Sm-^4 layering structure in an ori- 
entationally ordered phase. The shapes of V(r 12 ,Ui,u 2 ) 
for typical configurations are shown in Fig. 2. 

3. Simulation 

3. 1 Bulk phases 

The main purpose of this section is the determina- 
tion of the phases coexisting with the vapor phase. Here- 
after energies are measured in units of £1. The length-to- 
breadth ratio L/D is 2 throughout the paper. Although 
there are isotropic and solid phases but the Sm-^4 phase 
does not exist in the case that the hard SPC molecules 
with L/D — 2, 17 we expect that the additional orienta- 
tional interaction V may induce the Sm-^4 phase. The 
parameters e 2 and £3 are restricted in the region in which 
V is positive for any configuration of molecules. 

In performing the following simulations, a cutoff for in- 
termolecular interactions is introduced by means of shift- 
and-truncation. 22 The shifted-and-truncated potential is 
obtained by shifting the original potential by an amount 
S, followed by the truncation if the shifted potential is at- 
tractive and positive. The shift used in the present work 
is 5/ei = 0.01. The resulting cutoff distance is by defi- 
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Table I. The phases coexisting with a vapor phase in the Gibbs 
ensemble method. The letters I, N, and A denote the isotropic 
liquid, ncmatic, and Sm-/1 phases, respectively. 
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nition anisotropic, e.g., the cutoff distance in units of D 
is 4.075 in the end-to-end configuration and 3.120 in the 
side-by-side configuration for (£2, £3) = (0.6,-1.4). The 
maximum cutoff distance among the models that I have 
investigated is 5.097 for (£2, £3) = (1.2,0.0). 

The Gibbs ensemble Monte Carlo (GEMC) simula- 
tion 22 is used in this section since the GEMC simula- 
tion is appropriate to find a coexistence condition of two 
phases for the following reason. Two simulation boxes 
are used in the GEMC simulation. At a temperature suf- 
ficiently below the critical temperature at which the co- 
existence line of two phases ends, the spontaneous phase 
separation occurs in the GEMC simulation and the two 
simulation boxes are filled with these coexisting phases. 
Since the interface, whose effects on coexistence proper- 
ties are non-negligible, does not exist in GEMC simu- 
lation, quite large systems are not necessary for phase 
equilibrium calculations. Hence we can reduce the com- 
puter time utilizing the GEMC method. 

The details of the simulations are as follows. In all sim- 
ulation runs in this section, the total number of molecules 
is 512; the initial dimension of each simulation box is 
defined so that the reduced density of molecules is 0.1. 
Here the reduced density 77 is defined as 77 = N/(p c V) 
with the number of molecules N, the volume of the sim- 
ulation box V, and the hexagonal-close-packing density 
p c D 3 = 2/(V2 + V3L/D) of perfectly aligned hard SPCs. 
Each simulation box is maintained to be a cube through- 
out the simulation, and the periodic boundary condition 
is used. In order to make a uniform and isotropic ini- 
tial configuration, I performed a simulation for several 
thousands of MC steps at T = k^,T je\ = 2.0 (fee is the 
Boltzmann constant, T is the temperature, and T is the 
dimensionless temperature). The simulations were per- 
formed from higher temperature to lower temperature; 
the last configuration of a simulation at higher temper- 
ature was used as the initial configuration of a lower 
temperature simulation. In order to determine the upper 
and lower bounds of the transition temperature, heating 
processes are also simulated around the transition tem- 
perature, i.e., the last configuration at the lower tem- 
perature is used as an initial configuration of a higher 
temperature simulation. The sets of potential parame- 
ters investigated in this work are £2 = {0,0.6, 1.2} and 
£3 = {0,-0.5,-0.9,-1.2,-1.4}; £ 3 = -1.2 and -1.4 
are omitted at £2 = and 1.2 because V Q {r\2, ui, U2) 
takes a negative value depending on the configuration of 
molecules. 

The GEMC simulation results of the phases which co- 
exist with vapor phase are summarized in the Table I. 
These results show that £ 2 is necessary for the orienta- 
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Fig. 3. Phase diagram of LJ-SPC model with (£2, £3) = 
(0.6, —1.2) (upper panel) and (0.6, —1.4) (lower panel). Solid and 
open circles denote the reduced density observed in the cooling 
process and heating process, respectively. The vapor, isotropic 
liquid and Sm-yl phases are denoted by V, I and A, respectively. 
Data points of reduced densities rj obtained from four indepen- 
dent simulation runs are plotted at each temperature. 



tional ordering. The LJ-SPC model always exhibits the 
Sm-A phase if £ 2 is non-zero. The third parameter £ 3 is 
responsible whether the I or N phase is dominant. 

The LJ-SPC model with parameters (£2, £3) = 
(0.6,-1.2) and (0.6,-1.4), with which the model fluid 
undergoes a direct I- A phase transition, is of interest 
since the most of LT transitions observed in such LCs. 4 ~ 9 
The phase diagrams of the LJ-SPC model with these 
parameters are shown in Fig. 3. Each data point plot- 
ted in Fig. 3 is an average of r\ over 5 x 10 5 MC steps. 
The transition temperature Tja of the LJ-SPC model 
is, according to the Fig. 3, 0.597 < f IA < 0.610 for 
(e 2l s 3 ) = (0.6,-1.2) and 0.560 < f IA < 0.577 for 
(£2, £3) = (0.6,-1.4). For both sets of parameters, the 
Sm-^4 layer spacing at Tja is approximately 3 in units of 
D, i.e., approximately the molecular length L + D. 

3. 2 Smectic-A free standing film 

In this section, I determine a set of parameters for the 
existence of the Sm-^4 FSFs of LJ-SPC fluid, and investi- 
gate their properties. I also establish whether the LJ-SPC 
fluid corresponds to the LT materials. However, since the 
nucleation and growth of the dislocation loop are sup- 
pressed in a finite system, 23,24 we cannot expect to ob- 
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serve the LT without enormous amount of time. Thus I 
do not intend to observe the nucleation and growth of the 
dislocation loop, but focus only on the static properties 
characteristic to the LT materials. 

For simulations of the LC film in the NVT ensem- 
ble, I used an ordinary Metropolis method. The simu- 
lation box is a three-dimensional box with dimensions 
(L x , L y , L z ) = (15,15,50) in units of D with periodic 
boundary conditions in all directions. For the simula- 
tions in the present section, I used two kinds of initial 
state, which are prepared in the following ways: Sm-A 
initial state: slicing the Sm-^4 and vapor phase obtained 
from the GEMC simulation in §3.1 and sandwiching the 
Sm-j4 slice between the two vapor slices; uniform initial 
state: equilibrating the fluid at sufficiently high temper- 
ature (T = 2.0) for several thousands of MC steps until 
the fluid becomes uniform and isotropic. 

In the case that (£2, £3) = (0.6,-1.2) at the lower 
bound of the transition temperature T — 0.597 a sim- 
ulation from the Sm-^4 initial state shows that the initial 
Sm-^4 FSF melts into an isotropic film after 5 x 10 5 MC 
steps. From the uniform initial state at the same tem- 
perature, although the uniform isotropic fluid condenses 
into one or two films in equilibrium with the vapor phase, 
the condensed phase is still isotropic. The similar results 
are observed at T = 0.580 < Tia, hence Sm-^4 film is 
not stable in f > 0.580. At f = 0.570, the Sm-A FSF is 
stable. These results mean that the melting temperature 
of the Sm-^4 FSF is fairly lower than the bulk I- A phase 
transition temperature for (£2,23) = (0.6,-1.2). 

In the case (£2, £3) = (0.6, —1.4), from the Sm-^4 initial 
state, the Sm-^4 film does not melt at the lower bound 
of the phase transition temperature T = 0.560, unlike in 
the case £3 = —1.2. Furthermore, in the simulations at 
T = 0.560 from a uniform initial state, the LJ-SPC fluid 
spontaneously condenses into one or two Sm-^4 FSFs. 
One of a typical snapshot of the Sm-^4 FSF is shown in 
Fig. 4. The Sm-A FSFs with 3, 4, 5, and 6 layers are 
stable, but 1- and 2-layer Sm-^4 FSFs show strong fluc- 
tuation and seem to be unstable. The Sm-yl layer spacing 
is almost the same as that of the bulk Sm-yl phase. The 
Sm-^4 FSF prepared at T = 0.560 melts into isotropic 
film if it is heated up to T = 0.567. Since the hysteretic 
region of the Sm-^4 FSF is included in that of the bulk 
Sm-^4, we can expect that the I- A phase transition tem- 
perature of the Sm-^4 FSF is lower than or in the vicinity 
of the bulk Tia- Hence we can conclude that the LJ-SPC 
model with (£2, £3) = (0.6, —1.4) is similar to the normal 
Sm-^4 materials, which can produce the FSF in the vicin- 
ity of Tia, but the LJ-SPC model does not correspond to 
the LT materials, which exhibit the Sm-^4 FSF at higher 
temperature than the bulk Tia- In the rest of the present 
paper, I investigate the properties of the Sm-A FSF of 
the LJ-SPC model with (£ 2 ,£ 3 ) = (0.6,-1.4). 

In the case that the Sm-^4 initial state is used, the 
resulting condensed phase can be one or two films de- 
pending on the initial configuration and random-number 
sequence used in MC simulations. When the double-film 
state was obtained, I repeatedly performed the simula- 
tion run using different random-number sequences until 
a single-film state was obtained. The following numerical 
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Fig. 4. A snapshot of a 4-layer Sm-A FSF at T = T IA = 
0.560 with the potential parameters (£2,^3) = (0.6,-1.4). Each 
molecule is drawn as a segment of length 2D. 

results are obtained from single-film states and the data 
from double-film states are omitted. 

In order to check the SESO, I calculated a local density 
and local orientational order parameter. Here the local 
density p(z) is the number of molecules in [z, z + Az] 
divided with a volume L x L y Az. The local orientational 
parameter s(z) is defined as the maximum eigenvalue of 
an order parameter tensor Q(z) with components 

1 {Z) 1 
Q»v{z) = —^{Ui^Uiv) - -5^, fi,v = x,y,z (3) 

z i 

where the summation is performed over all the molecules 
within [z,z + Az] and N z denotes the number of 
molecules within [z, z + Az). In this paper, Az/D is 0.25. 
Figure 5 shows the local density and local orientational 
order parameter along the z axis at T = 0.556. Both 
the local density and local orientational order have max- 
imum peaks at innermost layers and gradually decrease 
to the free-surfaces, i.e., the free-surfaces reduce the Sm- 
A order in the LJ-SPC model. This is in contrast to the 
fact that the SESO is observed in experiments. 1, 2,9 

The disordering caused by the existence of free sur- 
faces can also be seen in the behavior of the averaged 
orientational order parameter s := J dz s(z)/L z . The N 
dependence of s is shown in Fig. 6. The data in this figure 
are obtained using ten thermal averages of s calculated 
from ten runs of 2 x 10 6 MC steps for each N, follow- 
ing an equilibration of 10 7 MC steps. The standard error 
of the mean is obtained from these ten averages with 
an assumption that the error distribution is a Gaussian. 
Figure 6 shows that s increases as the number of layers 
increases. This fact is another evidence for that the disor- 
dering effect of free surfaces strongly affects the thinner 
Sm-^4 FSFs and makes their order parameter small. 

The rest of this section is dedicated to the calculation 
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Fig. 5. (a)Local density p(z) and (b) local oricntational order 
s(z). The set of interaction parameters is (£2,£3) = (0.6,-1.4) 
and the temperature is T = 0.556. 
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Fig. 6. The number of molecules (TV) dependence of the averaged 
oricntational order parameter s at T = 0.556. The observed Sm- 
A FSFs consist of 3, 4, 5, and 6 layers in the case that N 6 
[470,530], [612,680], [740,820], and [872,930], respectively. 



of the film tension of the Sm-j4 FSF. The film tension: 

7= / dz[ Pn (z)-p t (z)}, (4) 
Jo 

is derived using the normal and tangential components 
of the pressure tensor 13, 14 defined respectively as: 



1 v-,„ dVij 



Ozij 
h3 J 



Pn (z) = p(z)f - —J2 
Pt{z) =p{z)f- t^5Z' i -,j 



(5) 



t,3 



xa-^- + Vio^- ) i (6) 



dxi 



dyi 



where (iy , ytj , z%j ) = Vj — r j , Vij is the potential energy 
between the ith and jth molecules, and V c is the volume 
of a slice [z, z+Az], i.e. V c = L x L y Az. The summation in 
eq. (6) is done if at least one of the ith and jth molecules 
is in the slice. 
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Fig. 7. The film tension 7 of Sm-A FSF with (e 2 ,£2) = 
(0.6, —1.4). The number of molecules is 612, with which the Sm- 
A FSF consists of 4 layers. 
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Fig. 8. The film tension 7 of Sm-A FSF with (£2, £2) = 
(0.6, —1.4) at T = 0.556. The simulation runs used in this figure 
are the same as those in Fig. 6. 



I measured first the temperature dependence of the 
film tension at some temperatures below the Tia- The 
data used in this calculation arc from MC runs per- 
formed in the same manner to the s calculation explained 
above (i.e., 10 7 MC steps for equilibration and ten runs 
of 2 x 10 6 MC steps for obtaining thermal averages). The 
temperature dependence of the average film tension and 
its standard error of the mean are shown in Fig. 7; the 
standard error of the mean is rather large as a result of 
the strong fluctuation of the film tension in simulation 
runs. Although it is difficult to state a definite conclu- 
sion because of the large standard error, the dependence 
of the film tension on T is very weak. 

From the calculation of the film tension with many 
initial configurations, it is found that the film tension 
strongly depends on the number of molecules N. The N 
dependence of the film tension is shown in Fig. 8. The 
set of simulation runs used in this figure is the same as 
that used in calculating s in Fig. 6. The Fig. 8 clearly 
shows that 7 decreases as N increases if the number of 
layers remains constant, and 7 increases discontinuously 
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at the point at which the number of layers increases. 

4. Concluding Remarks 

In this paper, I proposed a new LC model (LJ-SPC 
model) and showed that the model exhibits the bulk Sm- 
A phase and Sm-^4 FSF. The main features of the LJ- 
SPC model is the SPC-symmetric part and additional 
orientational part V a . The former is introduced to in- 
crease the stability of the Sm-^4 phase. The latter has the 
Maire-Saupe type interaction, which is often used in the 
molecular theoretical studies 25-27 to induce the orienta- 
tional order, and the term that governs the coupling be- 
tween translational and orientational degrees of freedom. 
In contrast to the models 25-27 neglecting the anisotropic 
repulsive part and the coupling between translational 
and orientational degrees of freedom, the LJ-SPC model 
exhibits the Sm-.A phase without resorting to the mean 
field approximation. 

From GEMC simulations in the bulk, the parameter 
£ 2 in eq. (2) is found to make a dominant contribution 
to the orientational order. The third parameter e 3 , which 
is set negative in this paper, contributes to diminish the 
orientational order in the bulk so that the N phase dis- 
appears if | £3 1 is large. The e 3 , on the other hand, plays a 
crucial role to stabilize the Sm-^4 FSF as shown in § 3.2. 
This is because the negative £3 term makes the side-by- 
side configuration stable and the end-to-end configura- 
tion relatively unstable, so that the perpendicular align- 
ment becomes to be preferred at a free surface; hence the 
Sm-^4 film becomes stable. 

It is clear that the LJ-SPC molecule is not appropri- 
ate as a model of the LT materials since the Sm-A FSF 
of the LJ-SPC is unstable above Tja- Furthermore, the 
experiments 9 show that the SESO is essential for the 
LT, but the free-surface of the LJ-SPC fluid reduces the 
Sm-^4 ordering. Hence we can conclude that the LJ-SPC 
model does not satisfy the condition to be a model of LT 
materials. 

The LJ-SPC model, however, bears a resemblance to 
the LT materials in the behavior of the film tension, 
which was precisely measured in ref. 7. The weak In- 
dependence of 7 below the Tja observed in ref. 7 is also 
observed in Fig. 7 in the present paper. Above the Tja 
the film tension observed in ref. 7 increases linearly with 
the temperature difference T — Tja at the constant num- 
ber of Sm-j4-layers, and discontinuously decreases if the 
layer-by-layer thinning occurs. Hence the Sm-^4 FSF of 
LT materials shows a sawlikc dependence of the film ten- 
sion on the temperature. The experiments in refs. 10 and 
9 showed that the melt of the innermost layer of the Sm- 
A FSF followed by the escape of the molecules into the 
meniscus occurs at a single layer-by-layer thinning transi- 
tion. It is natural to assume that the number of molecules 
in the FSF decreases monotonically as the temperature 
increases even in a temperature region between two suc- 
cessive layer- by- layer transitions. If that is the case, since 
the tension of the LJ-SPC FSF decreases as the number 
of particles decreases as shown in Fig. 8, we can assume 
that the sawlike dependence on the temperature of the 
tension of Sm-A FSF in ref. 7 is represented by the LJ- 
SPC FSF. 



The length-to-breadth ratio L/D, which is fixed at 2 in 
this paper, affects the stability of the Sm-^4 phase. The 
models with larger length-to-breadth ratio are important 
in considering the realistic Sm-yl LCs, since the real Sm- 
A LCs have much larger length-to breadth ratio (e.g. the 
4-cyano-4'-n-alkylbiphenyl (nCB) series exhibits Sm-^4 
phase only for n > 8). In general, the larger the length- 
to-breadth ratio is, the more stable the Sm-^4 phase be- 
comes. Hence we can expect that the bulk Sm-A phase 
and FSF will be stable for the LJ-SPC model not 

only with L/D = 2 but also with L/D > 2. In order 
to simulate the realistic Sm-j4, we have to adopt larger 
length-to-breadth ratios. Then the required values of £2 
and I £3 1 for the existence of the bulk Sm-^4 and Sm-^4 
FSF will be smaller than those of the LJ-SPC fluid with 
L/D = 2. That is, in the case of realistic Sm-^4 LCs, the 
higher order terms in the spherical harmonic expansion 
V will make a small contribution to the stability of the 
Sm-^4 phase. 

I add some comments about the absence of the SESO 
in the LJ-SPC fluid. In the presence of a V-A interface, 
the surface tension suppresses the layer fluctuation, 1-3 
which results from the Landau-Pierls instability. The re- 
sults in §IIIB, however, show that the effect of the sur- 
face tension is not enough to promote the SESO. The 
absence of the SESO in the present finite system may 
be attributed to the absence of a long-wavelength fluc- 
tuation, which gives dominant contribution to the total 
layer-fluctuation. 2,3 Since the SESO is induced by the 
suppression of the fluctuation at the surfaces, the under- 
estimate of the fluctuation may result in the omission of 
the SESO. There is a possibility that the SESO arises in 
the LJ-SPC FSF by including the long-wavelength fluc- 
tuation. The estimation of the effect of the long-wave 
length fluctuation on the SESO is left for future studies. 
In addition to such a hydrodynamic origin of the SESO, 
we can expect a molecular origin of the SESO. In order 
to promote the SESO, the authors in ref. 15 and 16 used 
a general spherical harmonic expansion of the attractive 
anisotropic interaction, and configured so that the terms 
of the interaction have the different ranges. Since the 
contribution of the terms favoring the Sm-^4 configura- 
tion to the free energy is proportional to the square of 
density gradient, these terms effectively deduce the total 
free energy if their ranges are sufficiently large to cover 
the interfacial region. Indeed, in ref. 15 and 16, setting 
the range of the terms favoring the Sm-^4 configuration 
largest, the authors showed by means of DFT that the 
SESO state has the minimum free energy. By the same 
reason, there is a possibility that the LJ-SPC fluid ex- 
hibits the SESO by introducing the concept of different 
interaction ranges. This requires further investigations. 

The negative film tension (i.e. interfacial free energy 
per unit area) observed in Fig. 8 requires further consid- 
eration. Since a negative interfacial free energy induces 
an expansion of the interface, the area of the interface 
tends to diverge. Therefore the negative film tension is 
inhibited to the thermodynamical stable state. However, 
in the present case, the area of the vapor and Sm-j4 in- 
terface cannot diverge without a distortion. Since the 
total free energy increases as a result of the distortion, 
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Fig. 9. (color online) The density profiles of 4-layer film with the 
number of molecules N = 470 (black solid curve) and 530 (red 
dashed curve). 



the distortion will be stopped at a state at which the 
free energy of the distortion and interfacial free energy 
balance each other. The distortion of the film can be ob- 
served in Fig. 9 which shows the 3-layer density profiles 
of N = 470 and 530; the peaks in the density profile of 
N = 530 film is broader than that of N = 470 film due 
to the film distortion. The averaged orientational order 
parameter s is also useful to estimate the distortion. For 
a fixed number of layers, s is an essentially decreasing 
function of the number of molecules as shown in Fig. 6. 
This fact indicates that the Sm-^4 FSF is distorted as a 
result of the negative interfacial free energy and there- 
fore the orientational order is reduced. Under conditions 
of the real experiment, however, the number of molecules 
consisting the FSF may decrease until the film tension 
becomes positive since the molecules can flow into the 
meniscus. Hence the distortion of the FSF may not be 
observed. 
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